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1 Introduction. 

The subject of the talk is the gluon structure function for DIS with nucleus. The gluon 
structure function is the most important physical observable that governs the physics 
at high energy (low Bjorken x) in the DIS. Dealing with nucleus we have to take into 
account the shadowing correction, which is the main point of interest in this talk. We 
show that the shadowing correction in the region of small x can be treated theoretically 
in QCD and can be reliably calculated using the information on the behaviour of the 
gluon structure function for the nucleon. We organize the presentation in the following 
way: first, we discuss the Glauber approach to the nucleus gluon structure function and 
answer the question what information on nucleon structure function we need to provide a 
reliable calculation using the Glauber formula; second, we briefly consider the corrections 
to Glauber approach that have been anticipated in QCD. It should be stressed that this is 
the first presentation of our results and the lack of space does not allow us to discuss the 
issue in details. This is why we are going to outline our strategy and to present the first 
estimates rather than to give the complete study of the problem which will be published 
elsewhere. 



2 Glauber approach in QCD . 

The idea how to write the Glauber formula in QCD has been first formulated in ref. 
M and was carefully developed by Mueller in ref. 0. It is easier to explain the idea 
considering the penetration of quark- ant iquark pair through the target. Indeed, during 
the time of passage through the target the transverse distance r t between quark and 
antiquark can vary by the amount Ar f oc R where E is the energy of the pair and R 
is the size of the target (see Fig.l). The quark transverse momentum (k t ) is k t oc ^- due 
to uncertainty principle. Therefore 




(i) 



holds if 



s > 



2mR 



(2) 



In terms of Bjorken x the above condition looks as follows: 
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(3) 



x < 



2mR 
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It means that the transverse distance between quark and antiquark is a good degree of 
freedom JI]] fl2] || . As has been shown by Mueller not only quark - antiquark pair can be 
considered in such way. The propagation of a gluon through the target can be treated in 
a similar way as the interaction of gluon - gluon pair with definite transverse separation 
r t with the target. The total cross section of the absorbtion of gluon(G*) with virtuality 
Q 2 and Bjorken x can be written in the form: 

a G * = (4) 
fdz f^p-j ^f(Q 2 ,r t ,x,z) 2 • {1 - exp[- a(r 2 ) S(b 2 t ) ] } • tff \Q\ r t , x, z) 

JO J ZTT J ZTT 

where i s the wave function of the virtual gluon with transverse polarization. As was 




Figure 1: The structure of the parton cascade in the Glauber formula. 
shown in ref. within leading log approximation of perturbative QCD (pQCD) we can 
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safely replace this function by t after integration over z in eq.(3). Finally the Glauber 
formula for the gluon structure function reads ( for N c = Nf = 3): 



where 

a GG = 3 -^^rUG(- 2 ,x) (6) 
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and S(b 2 ) is the profile function in impact parameter space for the interaction of the gluon- 
gluon pair with the target. We use for the calculation the Gaussian parameterization for 
5*, namely: 

S(tf) = ^ (7) 

where A is the number of the nucleons in a nucleus and R\ is the mean radius of a nucleus, 
which is equal to 

Rws is the size of the nucleus in the Wood - Saxon parameterization, which we chose 
Rws — tqAz with r = 1.3/m. Using the Gaussian parameterization for S we can take 
the integral over h t and get the answer: 



xG{Q\x) = 4 ^ ±J±f%{C + lnK G (x',r' t 2 ) + E 1 (k g (x' , rf) } (8) 

7T J x X J — 7t 7T T+ 

where C is the Euler constant and E\ is the exponential integral (see ref.0] 5.1.11) and 

K G (x',r' t 2 ) = 3a ° A Z r? x'G N (x', 4) (9) 



3 Theory status of the Glauber formula. 

In this section we would like to recall the main assumptions that have been made to get 
the Glauber formula: 

1. Energy (x) should be so high (small) to satisfy eqs.(2) and (3) and a s \n(l/x) ~ 1. 
The last condition means that we are doing the calculation in leading log(l/a;) approxi- 
mation of perturbative QCD (pQCD). 
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2. The GLAP || evolution equation holds in the region of small x. It means that 
a s \n{l/rf) ~ 1. One of the lessons that we have learned at this workshop is the fact 
that the GLAP equation is able to describe the HERA data quite well. 

3. Only the fastest partons ( GG pair) interacts with the target and there are no 
correlations between partons from different parton cascades (see Fig.l). 

4. There are no correlations between different nucleons in a nucleus. 

5. The average bt for GG pair - nucleon interaction is much smaller than Ra- 

We are going to discuss how well all the above assumptions work in the last section of 
the talk. 



4 Results. 

In our calculations we use the GRV parameterization |J for the nucleon gluon structure 
function. This parameterization describes the data quite well and it is suited for our 
purpose because (i) the initial virtuality for the GLAP equation is small and we can discuss 
the contribution of the large distances having some support in the experimental data; (ii) 
the parameterization uses the GLAP equation and the most essential contribution comes 
from the region where a s In Q 2 w 1 and a s ln(l/a;) ~ 1. 

4.1 Where the shadowing corrections are big. 

Fig. 2 shows the kinematic region of the deeply inelastic scattering. The curves are the 
solution of the equation kq — 1 for N (nucleon), Ca and Au. Above each of these curves 
the value of kg > 1 and the shadowing correction (SC) are big, below kq < 1 and the 
SC are rather small. 

4.2 What we are able to calculate in QCD. 

From the master equation (5) one can see that the large distances contribute to the value 
of the gluon structure function. Such contributions we are not able to calculate in pQCD 
and the value of the gluon structure function crucially depends on the hypothesis about 
nonperturbative behaviour of the gluon structure function that we have to assume to treat 
the large distances contribution. In pQCD we can safely calculate only the difference 
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Figure 2: Solution for k — 1. 
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xGa(x, Q 2 ) — xGa(x, Q 2 = Qq) where Ql is the initial virtuality. In Fig. 3 one can find 
the calculation for the ratio: 

R xG A (x,Q 2 ) - xG A (x,Q 2 = Q 2 ) 

1 A{xG N {x,Q 2 ) - xG N {x,Q 2 = QD) 1 } 

as function of x for Ca and Au (Q% = 0.25 GeV 2 ). 
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Figure 3: Ri as a function of ln{l/ xb) for Ca and Au. 
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4.3 Contribution of the large distances. 



As has been mentioned we are able to treat this problem only using some model for large 
distance behaviour of xG^. Fig.4 shows the ratio: 

2 A(xG N (x,Q 2 )) 1 ] 

for two models: 

1. xGn(x,Q 2 < Qq) = ^ xGn(x,Q 2 = Qq). This model takes into account the 

correct limit of the gluon structure function at small value of Q 2 , which follows from the 
gauge invariance of QCD. 

2. xGn{x,Q 2 < Qq) = xGn(x,Q 2 = Qq). In this model we assume that the scale 
for the behaviour xGn oc Q 2 is much smaller than Qq. We look at the difference in the 
value of i?2 as the estimates of possible errors that originated from our poor knowledge of 
long distance behaviour of the gluon structure function. The conclusion is that we cannot 
calculate the value of R2 even at x = 10~ 3 at Q 2 = lGeV 2 with better accuracy that 
20%, while at larger value of Q 2 (Q 2 ~ lOGeV 2 ) the accuracy is better ( about 5%). 



5 Correction to the Glauber formula. 

To abandon the main assumptions which have been made in the Glauber formula we have 
to develop a technique to include (i) the interactions of all partons ( not only the fastest 
one) with a nucleus; (ii) the parton interaction inside a nucleon and (iii) the nucleon 
correlation inside a nucleus. Such a technique has been suggested in ref . || and it is based 
on new evolution equation that takes into account the parton interaction inside the parton 
cascade as well as the parton interaction with the different nucleons. The lack of space 
does not allow us to discuss this problem in details but we want to point out that the 
Glauber formula shall be used as initial condition to the new evolution equation of ref. || . 

The numerical estimates |§ shows that the most essential contribution at least in 
HERA kinematic region is generated by the interaction of all partons with the target 
which corresponds to so called "fan" diagrams (see ref. 0) while the dynamic correlation 
(see refs.|| |I0| ) due to the interaction inside the parton cascade remains to be rather 
small. 

However, reliable estimates can be done only after extracting from HERA data the 
value of the scale for the SC for nucleon structure function. In our numerical estimates 
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which will be published elsewhere we follow the following strategy: we neglect the dynamic 
gluon correlations and iterate the master equation (5) several times. Since there is strong 
ordering in rapidity of parton in each parton cascade the "i — th" iteration means that we 
take into account the interaction with the target of all partons with the value of rapidity 
(y) larger than y^. It turns out that we have to make only two iterations for x > 10 -3 to 
get convergent result. 

6 Conclusions. 

We know the Glauber formula and the technique how to find the corrections to the Glauber 
formula in QCD. However we cannot provide reliable predictions for the gluon structure 
function for nucleus until we will get more data on low Q 2 and low x behaviour of the 
nucleon structure function. Unfortunately, we have to know not only the behaviour of the 
nucleon structure function but also extract from the experimental data the scale for the 
shadowing corrections to the value of the gluon structure function in the nucleon. We are 
going to check how the information from DIS with nucleus can reduce these uncertainties. 
At the moment we suggest to measure the ratio R\ which can be calculated with much 
better accuracy than the value of the gluon structure function. 
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